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ABSTRACT 

We present an analysis of the stellar populations of the gas-stripped outer disks of ten Virgo Clus- 
ter spiral galaxies, utilizing SparsePak integral field spectroscopy on the WIYN 3.5m telescope and 
GALEX UV photometry. The galaxies in our sample show evidence for being gas-stripped spiral 
galaxies, with star formation within a truncation radius, and a passive population beyond the trun- 
cation radius. We find that all of the galaxies with spatially truncated star formation have outer 
disk stellar populations consistent with star formation ending within the last 500 Myr. The synthe- 
sis of optical spectroscopy and GALEX observations demonstrate that star formation was relatively 
constant until the quenching time, after which the galaxies passively evolved. Large starbursts at 
the time of quenching are excluded for all galaxies, but there is evidence of a modest starburst in at 
least one galaxy. For approximately half of our galaxies, the timescales derived from our observations 
are consistent with galaxies being stripped in or near the cluster core, where simple ram-pressure 
estimates can explain the observed stripping. However, the other half of our sample galaxies were 
clearly stripped outside the cluster core. Such galaxies provide evidence that the intra-cluster medium 
is not static and smooth. For three of the most recently stripped galaxies, there are estimates for the 
stripping timescales from detailed gas stripping simulations. For all three of these galaxies, our strip- 
ping timescales agree with those from the gas stripping simulations, suggesting that star formation is 
quenched near the time of peak pressure. While the stripping of star-forming gas in the outer disk 
creates a passive population in our galaxies, there is still normal star formation in the center of our 
sample galaxies. It may be that Virgo is not massive enough to completely strip these spiral galaxies 
and, in a more dynamically active cluster or a cluster with a higher density ICM, such a process would 
lead to passive spirals and/or SOs. 

Subject headings: galaxies: clusters: individual (Virgo), galaxies: evolution, galaxies: interactions, 
galaxies: stellar content, intergalactic medium 



1. INTRODUCTION 

The morphology-density relationship (Dressier 1980) 
implies that clusters are creating SO galaxies at the ex- 
pense of spirals. The most salient difference between 
the two galaxy populations is the presence (spirals) or 
absence (SO's) of current strong star formation. Kauff- 
mann et al. (2004) find that the star formation - envi- 
ronment correlation is the strongest correlation of any 
galaxy property with environment, suggesting that envi- 
ronment plays a key role in the termination of star for- 
mation. While there is intense debate as to the specific 
environmental processes that cause the transformation, 
stellar population simulations (Shioya et al. 2004) have 
shown that simple quenching of star formation in a spi- 
ral galaxy can form SO-type spectra. However, not all 
SO galaxies are faded spirals; there is a significant frac- 
tion of SO's with larger bulges than spiral galaxies (e.g. 
Christlein & Zabludoff 2004), something that cannot be 
explained by simple disk fading. Ram pressure stripping 
is not the sole cause of the transformation of spirals to 
SO's, but it does appear to be a key part of the transfor- 
mation for many galaxies. 

A recently discovered population of cluster galaxies, 
"passive spirals" (Dressier et al. 1999), have noticeable 
spiral structure but little, if any, ongoing star formation. 

^ Current Address: Department of Astronomy, University of 
Massachusetts, 710 North Pleasant Street, Amherst, MA 01003- 
9305; hugh@astro.uniass.edu 



UV photometry (Moran et al. 2006) has shown that these 
galaxies have intermediate stellar population properties 
between spirals and SO's; passive spirals have no ongo- 
ing star formation, but have stopped forming stars only 
recently. These galaxies are consistent with an interme- 
diate phase between spirals and SOs. Poggianti et al. 
(2006) suggest that SO's may represent a heterogeneous 
population: "primordial passive galaxies" that formed 
at z > 2.5 and environmentally- formed SOs, which may 
evolve through a passive spiral phase. Shioya et al. (2004) 
demonstrate that such evolution from active spiral to 
passive spiral to SO is naturally explained by passive stel- 
lar evolution after a stripping event and also suggest that 
there are at least two paths to SO galaxies. 

In the nearby Virgo cluster, there exists a population 
of spiral galaxies with normal star formation inside some 
radius, but little or no star formation beyond that "trun- 
cation radius" (Koopmann & Kenney 2004). Many of 
these galaxies have relatively undisturbed stellar distri- 
butions, despite their severely truncated Ha disks. Such 
a disturbance in the gas with little or no apparent stellar 
disk disturbance strongly suggests a gas-gas interaction; 
if a gravitational interaction disturbed the star forma- 
tion, the stellar light should also be disturbed. The ap- 
pearance of these galaxies beyond the stripping radius 
(i.e. normal stellar disk with no current star formation) 
is strikingly similar to the passive spirals in higher red- 
shift clusters. By studying the stellar populations of our 
sample galaxies, we aim to understand when and where 
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in the cluster galaxies are stripped. Are these Virgo spi- 
rals with truncated star-forming disks just passive spirals 
that have not been fully stripped? 

Our sample of galaxies are drawn nearly entirely from 
the work of Koopmann, Kenney, & Young (2001), an 
atlas of _R-band and Ha-band images of spiral galaxies 
in Virgo. They find that ~ 50% of spiral galaxies in 
Virgo have truncated Ha disks (Koopmann & Kenney 
2004): relatively normal star formation in their central 
regions, with little or no star formation beyond a well- 
defined radius which we call the gas truncation radius. 
By studying the stellar populations of the outer disks 
of these galaxies through optical spectroscopy and UV 
imaging, we aim to understand where in clusters galaxies 
can be stripped and what effect gas stripping has on the 
star formation properties of spiral galaxies. The location 
in the cluster of stripping tells us how effective stripping 
is at different distances from the cluster center. Simple 
models of ISM-ICM stripping suggest that stripping can 
only occur in the cluster core, where the ICM densities 
are highest. Galaxies stripped far from the center of the 
cluster suggest a non-static, non-smooth IGM or extreme 
conditions outside the cluster core. 

In §2, we outline our observations and explain the ob- 
servational setup for optical spectroscopy and GALEX 
UV imaging. In §3, we discuss the relevant spectral fea- 
tures that we use to determine the stellar population 
properties of these galaxies. In §4, we describe the mod- 
els that we have used to interpret the stellar population 
properties of our observed sample. In §5, we detail our 
observations and results and discuss what the results im- 
ply about the quenching of star formation in the galaxies' 
outer disks. In §6, we discuss the implications of our indi- 
vidual galaxy results for the cluster as a whole, compare 
our results with simulations of neutral gas, and briefly 
discuss what our observations of these galaxies in Virgo 
may imply for galaxy evolution in clusters. Finally, in 
§7, we summarize our results. 

2. OBSERVATIONS 
2.1. Optical Spectroscopy 

Galaxies in the sample were observed the nights of 
March 28, 2003 - April 1, 2003 and March 18 - 23, 2004 
on the WIYN 3.5m telescope. For both runs, we used 
an identical spectrographic setup: the Bench Spectro- 
graphic camera with the 600@10.1 grating, tuned to pro- 
vide spectral coverage between 3900 A and 6800 A at 
a sampling of 1.4 A/pixel and a resolution of ^ 5.5 
A (FWHM). We used the SparsePak Formatted Field 
fiber instrument (Bershady et al 2004) in order to get 
large spatial coverage (80"x 80") of the Virgo galax- 
ies. While the spectral window extends from quite blue 
(~ 3900A) to fairly red (- 6800 A) the observing effi- 
ciency is not uniform over the entire range: fiber trans- 
mission causes significantly reduced throughput in the 
blue part of the spectrum and the grating sensitivity 
causes a drop in throughput in the red. 

Several factors were considered when deciding an ob- 
servational strategy for each of the disk galaxies in our 
sample. First, we must align the SparsePak array on the 
galaxy such that it's outside the region of significant Ha 
emission. This part of the galaxy, beyond the region of 
ongoing star formation, is interesting because it is a stel- 



lar clock telling us when star formation was interrupted 
in that region of the disk. Moreover, emission from star- 
forming regions would fill in the age-sensitive Balmer ab- 
sorption lines, leading to artificially old age estimates 
and complicating interpretations. Observationally, this 
makes it possible to observe the stellar populations in 
the disks of these galaxies by placing a majority of the 
SparsePak array outside the region of star formation (as 
shown in Figures 1-3). 

As an additional observational constraint, we observed 
as close to the Ha truncation radius as possible. This 
was done for both practical and scientific reasons. From 
an observational standpoint, the regions of the stellar 
disks are low surface brightness (often, at the Ha trun- 
cation radius, fiji > 21 mag/arcsec^), so it is necessary to 
observe in the highest surface brightness regions. From 
a scientific perspective, we are interested in when the 
galaxy was stripped of its star-forming gas. The inner 
regions should be the final regions to be stripped in an 
interaction with the ICM. Therefore, by studying the 
region of the galaxy close the stripping radius, we can 
determine how long ago the galaxy experienced strong 
pressure. If the galaxy is past the time of peak pressure 
stripping, this timescale will tell us how long ago peak 
pressure occurred. 

The SparsePak array was positioned such that the 
highest density of fibers is nearest to the truncation ra- 
dius. In most cases, the array was placed such that the 
central "diamond" of fibers was arranged along the ma- 
jor axis, just beyond the truncation radius. In some 
instances, however, the fiber array was placed parallel 
to the minor axis (i.e. NGG 4419, NGG 4569). In the 
case of NGG 4419, there is extended, low-level Ha emis- 
sion throughout the main disk, so we placed the array 
to attempt to avoid some of that emission. In the case 
of NGG 4569, we chose the location to ensure that we 
could effectively subtract sky light from the galaxy spec- 
tra. One potential complication associated with observ- 
ing along the minor axis is that bulge stars may con- 
taminate the light from the stellar disk. However, both 
of these galaxies have rather small bulges and in neither 
case do we use fibers placed on the center of the galaxy 
for our analysis. Therefore, we do not expect bulge light 
to significantly affect our results. 

The observations were reduced using the standard 
IRAF procedures, including zero correction and fiat field 
correction. Following basic image calibration, spectra 
were extracted using the dohydra task in IRAF, as well 
as wavelength calibrated and sky subtracted (using sev- 
eral SparsePak sky fibers) . The data were fiux calibrated 
using the spectrophotometric standard Feige 34. Finally, 
individual exposures are combined to create a set of 75 
spectra for each telescope pointing, with each spectrum 
being from one fiber (4'.' 7 or 360 pc^ in diameter). 

The spectra are then examined, fiber by fiber, for any 
residual emission that may be present in the regions of 
the galaxy observed. In most cases, it appears the Ha 
images (Koopmann, Kenney, & Young 2001) faithfully 
trace the emission and the spectra show little, if any, 
residual Ha emission beyond the truncation radius de- 
fined in the images. A comparison of two spectra from 

^ Here, and throughout this paper, we assume a distance to 
Virgo of 16 Mpc. 
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Fig. 1. — SparsePak Positions on R-band image (left) and Ha image (center). The composite spectrum from several summed fibers 
(indicated by the red circles on the images) is also shown (right). The radius given for each composite spectrum is the distance from the 
galaxy center to the center of the composite spectrum region. Shown here are images and spectra for NGC 4064, NGC 4388, NGC 4402, 
and NGC 4405. 
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Fig. 2.— Same as Figure 1, but for NGC 4419, NGC 4424, IC 3392, and NGC 4522. Note that the radial distance given for NGC 4419 
is distance along the minor axis from the galaxy's center. 



different regions of NGC 4522 shown in Figure 4 demon- 
strates tlie large differences between a spectrum taken 
beyond the truncation radius (top spectrum) and a spec- 
trum in the star- forming disk (bottom spectrum). The 
metal lines (Mgb and NaD) are largely unaffected, but 
the lower-order Balmer lines are clearly seen in emis- 
sion in the star-forming disk; in particular, H/3 is filled 
in by star-forming emission. While this effect decreases 
markedly for the higher order lines (Osterbrock 1989 re- 
ports emission line ratios of Ha/H/S = 2.86, Ha/H7 = 
6.09, Ha/H(5 =10.92 for a lO^K emission region), these 



lines still suffer from emission fill-in in the inner disk. 
While we went to some effort to exclude fibers with Ha 
emission from the composite spectra that we analyzed, 
it became apparent through the analysis that there was 
some minor residual Balmer emission in the outer disk 
in nearly all of the galaxies. The contamination is mi- 
nor, typically less than 0.25 A for H/? and much lower 
for H7 and US. Several fibers beyond the truncation ra- 
dius are combined for each galaxy to increase the signal 
to noise ratio; these composite spectra are used for all 
forthcoming spectral analysis. 
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Fig. 3. — Same as Figure 1, but for NGC 4569 and NGC 4580. Note that the radial distance given for NGC 4569 is distance along the 
minor axis from the galaxy's center. 



NGC 4522 

Major Axis 

r = 50" = 3.9 kpc 



il'"' ft 



%* 



Mgb 



NaD I 1 



H5 ay, 



HP 



Ha 



4000 4500 5000 5500 6000 6500 



NGC 4522 
Major Axis 
r = 20" = 1.5 kpc 



H8 Hy 



HP, Mgb , Nap 



Ha 



5000 5500 6000 

Wavelength (angstroms) 



Fig. 4. — Spectra from the disk of two regions of NGC 4522. The 
top spectrum is from Fiber 52, in the disk beyond the truncation 
radius (50" from the center and 10" from the H« truncation ra- 
dius). The bottom spectrum is from Fiber 53, closest to the center 
of the galaxy and still in the Ha emitting region. Note that there 
is no noticeable emission in the top spectrum and that Ha is seen 
in absorption. 
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For our data, which has intermediate spectral resolu- 
tion and intermediate signal-to-noise ratio, we have cho- 
sen to use the Lick/IDS spectral indices (Faber et al. 
1985). These indices are extracted from data by defining 
a spectral bandpass and two pseudo-continuum band- 
passes on either side of the feature of interest. 

We chose to extract Lick indices as opposed to directly 
measuring absorption lines because of the objectivity of 
the Lick system; while fitting absorption lines is appeal- 
ing for very strong features, the Lick system allows us 
to objectively extract spectral features without subjec- 
tively choosing which features of a line to fit. It should 
be noted that the indices that we extracted are not on 
the Lick/IDS system, since that would require smooth- 
ing our data to even lower spectral resolution. Instead, 
we leave our data at its native spectral resolution and 
just use the Lick index definitions (Worthey et al. 1994; 
Worthey & Ottaviani 1997) to compare our data with 
models. 

In practice, the spectra are analyzed using the the 
LECTOR task, which has been made publicly available 
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by A. Vazdekis'^. This Fortran code numerically inte- 
grates an input spectrum for several indices using Lick 
index definitions. Additionally, using S/N estimates of 
the spectrum, LECTOR measures the uncertainty in an 
index. The uncertainty in spectral indices that we quote 
throughout this work are the errors due to variations in 
the S/N and do not take into account any potential (un- 
known) calibration errors. 

2.2. GALEX UV Imaging 

For nine of our ten galaxies, we have also utilized 
GALEX imaging, from the 100s exposure All-Sky Imag- 
ing Survey (AIS) and from the ~ 1000s exposure Nearby 
Galaxy Survey (NGS) (Martin et al. 2005). More- 
over, we also have included, where possible, data from 
our medium-depth (1500 - 2500 s) Virgo GALEX sur- 
vey, which gives complete medium-depth coverage to the 
galaxies in the VIVA sample (Chung ct al 2008, in prep). 
The data include both far-UV (FUV; Aeff = 1516 A) and 
near-UV (NUV; Aoff — 2267 A) images and, in all cases, 
the galaxy was detected in the NUV. 

From this GALEX imaging, we extract fluxes in circu- 
lar apertures that match the size and location of the op- 
tical spectroscopy fibers. The appropriate apertures are 
then summed to determine the UV flux from the region 
of the composite spectrum. The FUV/NUV flux ratio 
for each galaxy is a further constraint on the stellar pop- 
ulation; because of how fast the FUV luminosity fades 
(see Figure 5), this measurement is extraordinarily sen- 
sitive to recent star formation for those galaxies where we 
are able to detect FUV emission. For the galaxies where 
FUV emission is undetected, the ratio of the FUV up- 
per limits (Morrissey et al. 2005) to the NUV flux gives 
an upper limit to the stellar population age. We have 
corrected all GALEX fluxes for Galactic absorption, us- 
ing the Schlegel et al. (1998) reddening maps and the 
extinction law of Li & Draine (2001). 

Due to uncertainties with pointing of ground-based 
spectroscopy we have gone to some care to make certain 
that we are sampling the same region with GALEX pho- 
tometry as we are with SparsePak spectroscopy. We have 
verified that the highest optical flux is present in the fiber 
centered at the nucleus of the galaxy in the SparsePak 
observations. Moreover, we have verified that strong Ha 
emission is not present in fibers that appear to be beyond 
the star- forming region in our images (Figures 1-3). In 
addition, we have verified that even sizable positional 
offsets do not significantly impact the color ratios. We 
have found that the flux ratios typically change by less 
than ±0.15 (0.05 in the log) for offsets up to -- 2.3" (the 
aperture radius). Moreover, in several cases, the colors 
are robust to significantly larger offsets. 

3. OPTICAL ABSORPTION LINE DIAGNOSTICS 

3.1. Age-Sensitive Indices 

The optical spectra that we use for our analysis are 
representative of the luminosity-weighted mean popula- 
tion for the region of the galaxy they are extracted from. 
This means that the observed spectrum is dominated by 
the younger, more luminous stellar population. While 
the extremely massive stars are the most luminous, stars 

^ See http://www.iac.es/galeria/vazdckis/vazdckis_softwarc.htrnl 



with spectral types O and B (with lifetimes of < 1 Myr 
and ~ 50 Myr, respectively) do not survive long. This 
means that, in a young population with no ongoing star 
formation, the stellar light is dominated by that from A 
stars for the majority of the first ~ 1 Gyr. 

The Balmer absorption lines have become the diagnos- 
tic of choice for measuring the age of the stellar popula- 
tion. For A-type stars or later, the strength of Balmer 
absorption lines are good indicators of stellar tempera- 
tures of stars at the main sequence turnoff and, there- 
fore, good indicators of the main sequence turn-off age. 
In addition, the Balmer lines are relatively insensitive to 
metallicity, allowing for a strong diagnostic not as con- 
fused by the age-metallicity degeneracy as some other 
indices. 

While we were careful to avoid most of the star- forming 
disk, in many cases, our absorption line measurements 
of H/3 are noticeably affected by nebular emission fill- 
in. Due to the superposition of absorption and emission 
features, it is difficult to quantify the degree of contami- 
nation from the emission regions. By comparing data to 
models, we are able to determine that the effect is small, 
in most cases < 0.5 A for 11(3, but it does noticeably 
affect the ages that one would infer from the H/3 index. 
Moreover, in all cases but one (NGC 4402), we have good 
measurements of the higher order Balmer lines (H7 and 
US). In fact, it is the ages inferred from the measure- 
ments of these lines that tell us that the H/3 line is con- 
taminated; in most cases, the age that we would infer 
from the H/3 line is older than the age inferred from the 
higher order lines. When comparing the observed spec- 
trum to a modeled spectrum, we do see a small amount 
of residual emission in the Ha feature which, assuming a 
Balmer decrement of Ha/H/3 = 2.86 (Osterbrock 1989), 
is enough to account for the age discrepancy between H/3 
and the higher order lines. Moreover, the age determined 
from H/3 is much more sensitive to errors in the measure- 
ments than the age estimates from H7 or H(5. Although 
the errors in H7 or US equivalent widths are higher than 
those in H/3 by as much as a factor of 3-4, the H7 and 
H(5 models cover ~ 5 times the range in their equivalent 
widths, compared to H/3. This makes the errors in the re- 
sulting ages somewhat smaller for the higher order lines. 
Additionally, because of the high ratio of Ha emission 
to H7 and US emission at (at typical HII region tem- 
peratures of ^ lO^K; Osterbrock 1989), there is much 
less Balmer emission in the higher order lines. There- 
fore, while we present the H/3 models for completeness, 
we use the average age concluded from the two higher 
order Balmer lines as our primary age diagnostic. 

3.2. Metallicity-Sensitive Index 

As our metallicity-sensitive index, we chose to use the 
hybrid index [MgFe]'. The index is a weighted average 
of the Mgb index and two Fe indices: 



[MgFe]' = VMgb * (0.72 * Fe5270 + 0.28 * Fe5335) (1) 

Thomas et al. (2003) specifically construct this index 
to be insensitive to variations in abundance ratios and 
a element enhancements for their sample of early type 
galaxies. While non-solar abundance ratios typically are 
not a concern for spiral galaxies, this index is still useful 
as an indicator of total metallicity. 



Stellar Populations of Spirals in Virgo 





TABLE 2 








GALEX Photometric Results 
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-16 
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4. STELLAR POPULATION MODELS 

Stellar population models describe how the light from 
a population of stars changes as they form and evolve. 
While models differ in their details, most share a com- 
mon, basic method of modeling a stellar population. 
Individual stellar evolution is modeled via isochrones 
from a stellar evolution code, a chosen IMF dictates 
the relative number of stars, the star formation history 
parametrizes how those stars formed as a function of time 
(i.e. SFR(t)), and empirical libraries or theoretical stellar 
atmospheres are used to transform the theoretical quan- 
tities of the isochrones (i.e. L and T^q) into observable 
quantities (i.e. broadband colors or spectra). 

4.L Starburst99 

For this work, we have chosen to use the Starburst99 
(SB99) stellar population synthesis code (Leitherer et 
al. 1999). Many stellar population synthesis codes use 
empirical stellar libraries to transform theoretical pa- 
rameters into observable quantities (i.e. Worthey 1994; 
Vazdekis 1999; Bruzual & Chariot 2003). A key short- 
coming of many of these spectral libraries is that they 
lack young, metal-poor stars (due to the dearth of these 
stars in the Solar Neighborhood). The outer spiral disks 
of this study contain young stellar populations. This 
is one of the primary reasons that we have selected the 
SB99 models to compare to the data: the most recent 
update of the SB99 models utilizes the synthetic stellar 
spectra of Martins et al. (2005), as opposed to empiri- 
cal spectral libraries. These models include metallicities 
from Z=0.001 (i.e. 1/20 solar) to Z=0.04 (i.e. - twice 
solar), allowing more complete metallicity coverage. This 
means that we avoid the problems of the incomplete li- 
braries and are able to define a uniform grid of models, 
from which we can extract ages and metallicities. 

4.2. Transformation to Observables 

As its output, SB99 produces spectra created by sum- 
ming several theoretical spectra of Martins et al. (2005). 
These spectra are then smoothed to the instrumental res- 
olution of the observed spectra. While the resolution of 
the spectra will not affect the equivalent width measure- 
ments, the Lick index measurements can be affected by 
spectral resolution mis-match since some line flux can be 
smoothed into the pseudo-continuum bandpasses. While 
this is not ideal as some of the line flux is shifted out 
of the line passband, as long as there is consistency be- 
tween the observations and models, our conclusions will 



be robust. Unlike many stellar population analyses of 
earlier-type galaxies, we do not correct our galaxies for 
stellar velocity broadening of absorption lines. Since we 
are studying the disks of spirals, the velocity dispersion is 
very low. In the case where the effect should be most ex- 
treme (NGC 4388, a massive edge-on galaxy), we find no 
meaningful difference in results from the corrected and 
uncorrected spectra. 

SB99 additionally generates lower resolution SEDs 
(Leitherer et al. 1999) that extend far into the UV. These 
SEDs are used, together with the detector bandpasses for 
GALEX to create a grid of colors that depends critically 
on age. Specifically, because the FUV fiux changes so 
quickly at young ages (Figure 5), the FUV to NUV flux 
ratio is a powerful indicator of population age. 

4.3. Modeled Star Formation History 

Much of the stellar population synthesis work in the lit- 
erature relies on the "Simple Stellar Population" (SSP) 
assumption; that is, the assumption that all stars formed 
at one time and then passively evolved. While this may 
be a good approximation of star clusters, it is clearly not 
a good description of the star formation history in a spiral 
galaxy; one of the fundamental characteristics of a spiral 
is ongoing star formation. Therefore, we have chosen to 
model the stars in the outer disks of these Virgo galaxies 
as quenched star formation histories. That is, we assume 
that there was constant star formation from a formation 
time (set to be 12 Gyr ago), followed by rapid quenching 
of star star formation. We define the time since star for- 
mation ceased as the quenching time (tq). Note that tq is 
different from the "age" of the stellar population, which 
is the value usually quoted when comparing observations 
to models. The luminosity-weighted mean age of our 
quenched stellar population is significantly larger than 

iq. 

Practically, the quenched star formation models are 
constructed by summing spectra from several bursts. 
Constant star formation is approximated as a star form- 
ing burst every 50 Myr from the population's formation 
time (12 Gyr) to the time of quenching (tq). In addition, 
we have considered the effect of adding a small starburst 
at the time of quenching, in which 2% of the stellar mass 
is formed in a burst at the time of quenching (equivalent 
to a ~ 5 times increase in SFR for 50 Myr). 

The composite theoretical spectra are analyzed in the 
same manner as the observed spectra. LECTOR is used 
to measure the absorption line values of each model, and 
the absorption line strengths extracted from these spec- 
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Fig. 5. — Evolution of the spectrum from SB99 quenched star 
formation models with constant star formation followed by instan- 
taneous quenching. Shown are models with quenching times of 50 
Myr, 250 Myr, 500 Myr and 1 Gyr ago. Also indicated on the 
figure are the FUV and NUV bandpasses for GALEX, along with 
the Johnson broadband V bandpass. 

tra allow us to form a line index grid (Figures 6 - 15) from 
which we can extract a value of tq for any given galaxy by 
determining the line of constant age that passes through 
the absorption line measurement. The errors in tq are 
then calculated by determining the line of constant age 
that pass through the absorption line index ± the error 
in the index. 

5. THE GALAXIES 

When wc introduce three groups of galaxies in the 
following section, we will quote the luminosities and 
galaxy morphological types (Table 3), in addition to 
discussing the star-forming properties. The luminosi- 
ties (in terms of L,) refer to the K band luminosities 
of each galaxy (as determined by 2MASS and assum- 
ing a distance of 16 Mpc). While we list the mor- 
phological types, it has been pointed out that these 
types do not work well in the Virgo cluster (Koop- 
mann & Kenney 1998). In particular, nearly all of the 
galaxies that we have studied have truncated Ha disks. 
The star formation properties are drawn from Koop- 
mann & Kenney (2004), which discusses the star form- 
ing properties for 55 Virgo Cluster spirals. Nearly all 
of our galaxies fall into their "Truncated/Normal" or 
"Truncated/Compact" categories. Galaxies classified as 
"Truncated/Normal" have no star formation beyond a 
"truncation radius" (varies by galaxy between 0.3-R25 - 
O.8/J25), but normal star formation interior to that ra- 
dius. Galaxies classified as "Truncated/Compact" also 
have an outer disk free from star formation, but have 
less regular inner disks, with non-axisymmetric star for- 
mation confined to the inner '^ 1 kpc. We also report the 
HI deficiency (Def///), as defined by Giovanelli & Haynes 
(1983) and as measured by Ghung et al. (2008). This 
quantity is the logarithmic ratio of the HI mass expected 
for an isolated spiral galaxy to the HI mass measured in 
each spiral galaxy. A Defni—O indicates an HI- normal 
galaxy, Defni — ^ indicates a galaxy that has 10 times less 
HI than expected, and so on. In Table 3, we also quote 
the projected distance of the galaxies from M87 (dsr), 
both in terms of measured angular separation and com- 
puted linear distance (assuming a distance of 16 Mpc). 

For all of the galaxies in our sample, we present a 
timescale that represents the time since quenching, as- 
suming a quenched star formation model- tq. The value 



of tq quoted is the average of the values determined from 
H5A and H7A measurements; the age determined from 
H/3 measurements tends to be unreliable due to some 
nebular emission fill-in. 

In Figures 6 - 15, we present the optical spectrum 
results for our sample of Virgo galaxies. Each figure 
presents the Balmer line strengths (H(5A, H7A, and H/3) 
plotted against the [MgFe] ' strength (these line strengths 
are also listed in Table 4 and the full spectra can be seen 
in Figures 1-3). We compare these line strengths to 
Starburst99 models for two different star formation his- 
tories. The top row of the figures shows an evolutionary 
track for a SSP model with solar metallicity, with some 
sample ages (in Gyr) marked along the line. While the 
star formation history of spiral galaxies is likely markedly 
different from an SSP, these models are shown here for 
comparison. Models are displayed in the middle rows 
for galaxies with a constant star formation rate to the 
quenching time, then no star formation for a range of 
metallicities: Z=0.008, Z=0.02 (solar), and Z=0.04. Fi- 
nally, in the bottom row, we show the region of the spec- 
tra with Balmer line absorption. In each case, we de- 
termine the time since star formation was quenched by 
comparing the measured line strengths to SB99 models 
(also see Tables 5 & 6). 

5.1. Truncated "Normal" Galaxies 

A majority of the galaxies in our sample have normal 
stellar disks and Ha disks that are truncated at ^ 0.3i?25 
- 0.8i?25. The stellar disks, as traced by the R band light 
are generally well-fit by elliptical isophotes and do not 
show any apparent stellar asymmetries. 

5.1.1. Old Outer Stellar Disks 

A subset of the Truncated "Normal" Galaxies are those 
with old outer stellar disks. The Ha morphology for 
NGG 4405 (Figure 6), IG 3392 (Figure 7), and NGG 4580 
(Figure 8) are virtually identical; all three galaxies have 
Ha disks that are small compared to their optical disks 
and dramatically truncated at a small fraction of their 
optical radius. NGG 4419 (Figure 9) also has a truncated 
Ha disk, but the existing Ha emission is much more ane- 
mic than the other three galaxies. All four of the galaxies 
have undisturbed outer stellar disks. The outer disk of 
NGG 4580 is somewhat unusual because, despite a lack 
of ongoing star formation, there are stellar spiral arms in 
the outer disk. 

For all four galaxies, the optical spectrum is relatively 
flat, but shows a significant break in the blue region, 
leading to an overall "red" color. Comparison of the 
spectra to models tells us that the star formation in all 
of these galaxies was quenched ~ 450 Myr ago or longer. 
None show any evidence for extraplanar HI or asymmet- 
ric star formation and all are highly HI deficient, with, at 
most, 10% of their "expected" HI content (Ghung et al. 
2008; Giovanelli & Haynes 1983). All of these galaxies 
have been observed with GALEX and none of them have 
significant FUV emission in the outer disk, as would be 
expected for regions of galaxies with no star formation 
within the last - 400 Myr. 

5.1.2. Young Outer Stellar Disks 

A second subset of those galaxies with undisturbed 
stellar disks are those that have young outer stellar disks. 
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TABLE 3 
Properties of Sample Galaxies 



Galaxy 1 


dfi7 (deg) 


^87 (Mpc) 


Vr (km/s) 


Type 


SF Classt L (L,,k) 


in 


Dcfff/t 


NGC 4064] 


8.8 


2.45 


913 


SBc(s) 


T/C 


0.6 


75 


1.75 


NGC 4388 


1.3 


0.35 


2524 


Sb 


T/N 


0.9 


82 


1.18 


NGC 4402 


1.4 


0.40 


239 


Sb 


T/N 


0.6 


80 


0.74 


NGC 4405 


4.0 


1.11 


1747 


Sc/SO 


T/N [s] 


0.3 


53 


0.98 


NGC 4419 


2.8 


0.78 


-196 


Sa 


T/A 


0.5 


82 


1.23 


NGC 4424 


3.1 


0.87 


440 


Sa pec 


T/C 


0.3 


63 


0.75 


IC 3392 


2.7 


0.75 


1683 


Sc/Sa 


T/N 


s 




0.3 


70 


1.16 


NGC 4522 


3.3 


0.92 


2328 


Sc/Sb 


T/N 


s 




0.5 


78 


0.86 


NGC 4569 


1.7 


0.47 


-220 


Sab(s)I-II 


T/N 


s 




3.4 


66 


1.52 


NGC 4580 


7.2 


2.0 


1033 


Sc/Sa 


T/N 


s 




0.5 


42 


1.44 



tStar Formation Class, as defined by Koopmann & Kenney (2004): (A) anemic, (T/A) 
truncated/anemic, (T/C) truncated/compact, (T/N) truncated/normal, (T/N [s]) trun- 
cated/normal (severe) 

■tHI Deficiency from Chung et al. (2008), as defined by Giovanelli & Haynes (1983). 



TABLE 4 
Absorption Line Strengths for Composite Spectra 



Galaxy 



H<5A (A) 



H7A (A) H/3 (A) Mgb (A) Fe5270 (A) Fe5335 (A) [MgFc]' (A) 



NGC 4064 
NGC 4388 
NGC 4402 
NGC 4405 
NGC 4419 
NGC 4424 
IC 3392 
NGC 4522 
NGC 4569 
NGC 4580 



2.15 ±0.68 
3.39 ±0.98 

1.97 ±0.80 
1.66 ±0.65 
3.20 ±0.71 
1.15 ±0.80 
5.65 ±0.64 
2.35 ±0.68 
-0.26 ±0.86 



-0.66 ±0.51 

1.62 ±0.76 

1.28 ±0.77 

-1.16 ±0.64 

-2.94 ±0.48 

0.88 ± 0.54 

-1.43 ±0.59 

4.18 ±0.51 

-1.10 ±0.50 

-3.39 ±0.65 



2.84 ± 0.24 
3.78 ±0.40 
2.57 ±0.39 
3.11 ±0.31 
1.97 ±0.21 
3.24 ±0.27 
2.80 ±0.28 
4.33 ±0.28 
2.64 ± 0.23 
2.48 ±0.29 



2.22 ±0.23 

2.33 ±0.41 
1.49 ±0.38 
2.49 ± 0.31 
2.94 ±0.19 
2.08 ±0.27 

2.34 ±0.27 
1.61 ±0.30 
2.80 ±0.23 

3.23 ± 0.28 



2.45 ±0.25 
2.15 ±0.44 
1.96 ±0.40 
2.26 ±0.33 
2.45 ±0.21 
2.15 ±0.29 
2.47 ±0.29 
1.78 ±0.32 
2.62 ±0.24 
2.83 ±0.30 



2.30 ±0.28 
1.37 ±0.50 
1.11 ±0.46 
1.90 ±0.37 

2.30 ±0.23 
1.85 ±0.32 

2.31 ±0.32 
1.58 ±0.37 
2.29 ±0.27 
2.53 ±0.33 



2.31 ±0.15 
2.12 ±0.26 
1.60 ±0.25 

2.32 ±0.20 

2.66 ±0.13 
2.07 ±0.18 
2.38 ±0.18 

1.67 ±0.20 
2.66 ±0.15 
2.97 ±0.18 



These include NGC 4388 (Figure 10), NGC 4522 (Fig- 
ure 11), and NGC 4569 (Figure 12). These galaxies aU 
have disturbed HI spatial distributions (Oosterloo & van 
Gorkom 2005; Kenney et al. 2004, 2008). The star forma- 
tion in all of these disks is highly truncated. NGC 4388 
has a small, relatively symmetric Ha disk, although there 
is a known Ha (Yoshida et al. 2002) and HI (Oosterloo 
& van Gorkom 2005) tail extending to the northeast. 
The star formation in the disk of NGC 4522 is strongly 
truncated, although there is significant extraplanar Ha 
and HI emission to the northwest of this galaxy (Kenney 
et al. 2004). NGC 4569 also has strongly truncated Ha 
emission, but has an HI and Ha arm that extends out 
from the NE edge of the truncated gas disk (see the Ha 
image in Figure 3) (VoUmer et al. 2004; Kenney et al. 
2008). 

The spectra of the outer disks of these galaxies are ei- 
ther moderately (NGC 4388, NGC 4569) or extremely 
(NGC 4522) blue, indicating young stellar populations. 
All of the galaxies show strong Balmer line absorption in 
the outer disk, with NGC 4522 having extremely strong 
absorption lines. In addition, all of these galaxies have 
moderately weak metal lines when compared to galax- 
ies with older stellar populations; this is expected in a 
population with a higher fraction of A- type stars. In 
addition, all of these galaxies have significant GALEX 
emission, both in the FUV and NUV. NGC 4522 is the 
extreme case once again, as the only galaxy outer disk 
with stronger FUV emission than NUV emission. By 
combining UV data and optical spectroscopy, we con- 



clude that star formation persisted in the stripped disk 
of NGC 4522 up to 50-100 Myr ago (Growl & Kenney 
2006). 

Star formation extends out to ^ O.7R25 in the disk in 
NGC 4402 (Figure 13), so the stellar surface brightness is 
relatively low beyond the Ha truncation radius. There- 
fore, the stellar population results are somewhat incon- 
clusive. The H/3 index seems to indicate a older popula- 
tion, but other evidence evidence is in conflict with this, 
perhaps because the H/3 line is partially filled in with 
emission. First, the GALEX photometry suggest that 
there is a young stellar population in the outer disk of 
this galaxy. Secondly, the H7 line is much stronger than 
would be expected for an old population (Figure 13). In 
addition, there is abundant evidence from optical imag- 
ing and radio mapping that this galaxy is currently un- 
dergoing stripping (Growl et al. 2005; VoUmer et al. 2007; 
Murphy et al. 2008). Unfortunately, we are unable to 
accurately measure the H(5 line in the outer disk of this 
galaxy because the surface brightness of the stellar disk 
at this location is too low. 

5.2. Galaxies with Complex Interaction Histories 

Finally, there are two galaxies in the sample that ap- 
pear to not have been simply stripped. In the case of 
NGC 4064, Cortes et al. (2006) note that outer isophotes 
appear undisturbed, but that there is intense circum- 
nucular star formation in a bar-like feature and dis- 
turbed central dust morphology. Cortes et al. (2006) 
interpret the broadband morphology. Ha morphology. 



10 



Crowl & Kenncy 



NGC 4405 



12-. 



H6A 



\ 3 



10 -5 



# 



0.5 ;; 



0.05 ^' 



12^ 



0.05^' 



HyA 


12 

\ 


Hp; 


- 


1 

'^3 


: 




^ 




-._ 0.5 


0.05 ~ 


-_ 0.5; 

_y : 



5 10 -10-5 5 1017345 

H6A H7A H(S 




4050 4150 

Observed Wavelength (Ang.) 



4300 4400 

Observed Wavelength (Ang.) 



4860 4900 

Observed Wavelength (Ang.) 



Fig. 6. — Absorption line indices measured for Hi5 (left column), H7 (middle column), and H/3 (right column). The top row shows 
comparison with a solar metallicity SSP, which indicates the luminosity-weighted mean age of the population. The middle row shows 
comparison with a quenched model, with tq indicated by vertical solid lines corresponding to 0.05, 0.5, 1.0 and 12.0 Gyr. The three broken 
lines correspond to metallicitics of Z=0.008 (short dash), Z=0.02 (i.e. solar; dot-dash), and Z=0.04 (long dash). The bottom row shows the 
region of the spectrum from which the Balmer lines were extracted with the index bandpass, the continua bandpasses, and the continuum 
line indicated. 



and disturbed gas distribution as a signature that this 
galaxy experieneed a gravitational interaction, in addi- 
tion to possibly being stripped of its gas through an in- 
teraction with the ICM. NGC 4424 also has compact 
central star formation in a bar-like feature. In con- 
trast to NGC 4064, the broadband optical appearance 
of NGC 4424 is also peculiar, with shell-like features 
and banana-shaped isophotes, suggesting a major grav- 
itational disturbance or a merger (Kenney et al. 1996). 
Additionally, the central 2-3 kpc show disturbed dust 
lanes and infalling molecular gas (Cortes et al. 2006). 

The spectrum of NGC 4064 shows characteristics of an 
intermediate-aged population (Figure 14). The spectrum 
is fairly red and shows strong line blanketing in the blue 
part of the spectrum. The metal lines (Mgb and NaD) 



are strong, indicative of a population dominated by giant 
stars. If we consider the quenching time as extracted 
from the US and II7 lines, we find tq ~ 425 ± 75 Myr. 
This age strongly suggests that NGC 4064 was stripped 
in the cluster outskirts and not in the core. See further 
discussion in Section 6.3. 

The stellar population in NGC 4424 is one of the 
younger in our sample (Figure 15). This fact is obvi- 
ous by looking at the composite spectrum, which has a 
moderately blue color, sloping downward to the red part 
of the spectrum. Additionally, the Balmer features are 
quite strong, with all three major lines easily distinguish- 
able. When comparing these data to models, we obtain 
a quenching time of tq = 275 ± 75 Myr for the higher 
order line indices. From the H/3 measurements, we infer 
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Fig. 7.— Same as Figure 6, but for IC 3392. 
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that there is some modest residual ionized gas emission 
in the outer disk, but H(5 and H7 give a consistent age, 
and the effect appears to be minor. 

GALEX observations of NGC 4424 indicate a mod- 
erately young disk and agree with the results from the 
stellar spectroscopy. The relatively blue GALEX colors 
indicate that star formation in this galaxy was quenched 
~ 300 Myr ago. The timescales from UV imaging and 
optical spectroscopy only agree for a quenched star for- 
mation history; therefore, these data are consistent with 
constant star formation, followed by an instantaneous 
quenching. 

6. DISCUSSION 

6.1. Ages and Star Formation Histories from GALEX 
Observations 

Because of the extraordinary sensitivity of GALEX UV 
observations to the presence of young populations, we 
can potentially distinguish between different star forma- 



TABLE 5 

Time Since End of Star Formation of Quenched 

Spirals as Determined for Individual Balmer 

Indices 



Galaxy 


tq,HiA(Gyr) 


tq,H7A(Gyr) 


tq,H/3(Gyr) 


NGC 4064 


0.4 ±0.1 


0.45 ±0.1 


1.5 ±0.9 


NGC 4388 


0.25 ±0.1 


0.2 ±0.1 


0.25 ± 0.2 


NGC 4402 




0.2 ±0.15 


12 


NGC 4405 


0.4 ±0.15 


0.5 ±0.15 


0.6 ±0.2 


NGC 4419 


0.35 ±0.1 


0.65 ±0.1 


> 12 


NGC 4424 


0.275 ± 0.075 


0.275 ± 0.075 


0.55 ±0.2 


IC 3392 


0.5 ±0.1 


0.5 ±0.1 


12+'' 


NGC 4522 


<0.05 


<0.05 


<0.05 


NGC 4569 


0.3 ±0.05 


0.3 ±0.05 


1 2+3 
-'-■^-0.5 


NGC 4580 


0.45 ±0.1 


0.5 ±0.1 


10+3 

-'-■'^-0.4 



tion histories using our GALEX observations. In partic- 
ular, we can extract information about a galaxy's star 
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TABLE 6 

Prcjected Distance from M87 

AND Derived Stellar Population 

Quenching Timescale for the 

Stripped Spirals 

Galaxy ds? (Mpc) iq (Myr) 



NGC 4064 


2.45 


425 ± 75 


NGC 4388 


0.35 


225 ± 100 


NGC 4402 


0.40 


< 200 


NGC 4405 


1.11 


450 ± 150 


NGC 4419 


0.78 


500 ± 150 


NGC 4424 


0.87 


275 ± 75 


IC 3392 


0.75 


500 ± 100 


NGC 4522 


0.92 


100 ± 50 


NGC 4569 


0.47 


300 ± 50 


NGC 4580 


2.01 


475 ± 100 



formation history by combining analyses of UV emis- 
sion, which traces a very young stellar population and 



stellar absorption lines, which give us information about 
star formation over a longer timescale. There is mod- 
erately good agreement between the ages determined 
by spectral analysis of the outer disk and the GALEX 
UV observations (Figure 16) when assuming a quench- 
ing model. There is a clear sequence from NGC 4522, 
which has strong Balmer lines (young ages) as deter- 
mined by spectral measurements and a very hard UV 
color, to NGC 4580 and IC 3392, which have compar- 
atively weak Balmer lines (older ages) and significantly 
softer UV colors. While there is a clear sequence, the 
galaxy measurements are clearly offset from the solar- 
metallicity models (top panel of Figure 16). As can be 
seen when comparing to the Z=0.008 models (bottom 
panel of Figure 16), this is apparently a metallicity effect, 
showing that there is some degeneracy between star for- 
mation histories (i.e. quenching, quenching plus a burst, 
SSP) and mctallicities. Metallicity estimates from the 
[MgFe]' index indicate that most of the galaxies have a 
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luminosity-weighted mean metallicity that is intermedi- 
ate between Z^^O.OOS and Z=0.02, suggesting that all of 
the galaxies are consistent with simple quenching of star 
formation, with the exception of NGC 4522 (See Growl 
& Kenney 2006). 

As discussed in Growl & Kenney (2006), the outer disk 
of NGC 4522 is most consistent with a small starburst 
at the end of its outer-disk star formation, such that 
2% of the total stellar mass is formed in a very short 
interval at the time of quenching. As shown in Figure 
11, the absorption line indices can not be explained with 
a quenched star formation history; the Balmer lines in 
the outer disk are too strong. It is through GALEX 
observations (shown in Figure 16 and in Growl & Kenney 
2006) that we are able to obtain a stripping timescale 
without making an assumption about the star formation 
history. 



6.2. Comparison of Stellar Population Ages to Gas 
Stripping Simulations 

6.2.1. Timescale Comparisons 

For three of the galaxies in our sample, we can compare 
our timescale estimates to estimates from a completely 
independent analysis. While we observe and model the 
stellar properties of the stars left behind after ISM-IGM 
stripping, it is also possible to measure and model the 
gas that has been stripped from these galaxies. For most 
galaxies in our sample, the lack of high column density 
gas near the disk of the galaxy means it is difficult to 
constrain the interaction parameters from gas modeling. 
However, for three galaxies, NGG 4522, NGG 4388, and 
NGG 4569, we are able to compare our results to those 
from gas stripping models. 

We find a very young star formation quenching age for 
NGG 4522, tq = 100 ± 50 Myr. This implies that star 
formation has ended very recently and, indeed, we do see 
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Fig. 10.— Same as Figure 14, but for NGC 4388. 



large amounts of cxtraplanar gas near the disk midplane 
indicating that stripping is ongoing. A comparison of HI 
distribution and kinematics with ram-pressure stripping 
simulations of this galaxy (VoUmer et al. 2006) imply 
that peak pressure was ^ 50 Myr ago, in excellent agree- 
ment with our results. 

For NGC 4388, the stellar population analysis indi- 
cates that this galaxy was stripped longer ago, tq — 225 
Myr. A comparison of the HI distribution and kine- 
matics with a ram pressure simulation of this galaxy 
(VoUmer & Huchtmeier 2003) originally suggested that 
it was stripped ^ 120 Myr ago. This is in reasonable 
agreement with our data, but shows some difference be- 
tween the HI simulations and the stellar population ages. 
However, recent discoveries of more HI much further from 
the galactic plane (Oosterloo & van Gorkom 2005) im- 
ply that the stripping must have happened longer ago. 
Specifically, considering how far the HI is from the galac- 
tic plane, the time since peak pressure should be revised 
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upward to ^ 200-250 Myr (B. VoUmer, private commu- 
nication). While these models are based on much less 
complete HI data than the models for NGC 4522 and 
NGC 4569, the stripping ages are still in agreement with 
our estimates from the stellar population analysis. 

For NGC 4569, we find the oldest stripping age of the 
three: tq ~ 300 Myr. This is precisely the age that 
VoUmer et al. (2004) derive from simulations of their HI 
observations of this galaxy. In fact, VoUmer et al. (2004) 
conclude that the HI features can only be explained if 
this galaxy is significantly past the peak pressure phase 
and much of the observed HI is falling back onto the disk. 
Additionally, Boselli et al. (2006) model the broadband 
color profile of NGC 4569 and find that the time since 
interaction is constrained to be less than 400 Myr ago. 

6.2.2. Implications of Timescale Comparisons 

The agreement of our results with those from stripping 
simulations (VoUmer & Huchtmeier 2003; VoUmer et al. 
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2004, 2006) is remarkable for several reasons. First, the 
agreement of our estimates for quenching timescales 
with their stripping timescales is encouraging. While 
both our stellar population analysis and HI morphol- 
ogy and kinematic analysis are built on a large set of 
assumptions, the assumptions are generally different for 
the two methods. While the stellar population analysis 
relies on assumptions about star formation history and 
stellar evolution models, the HI simulations rely on how 
the ICM interacts with the ISM, how gas is transported 
away from the disk and the ram pressure profile that a 
galaxy experiences while interacting with the ICM. In the 
case of these three galaxies, the agreement between the 
timescales from HI models and the timescales from our 
stellar population analysis suggests that our timescales 
are probably correct. 

Secondly, the agreement between the ages derived from 
HI data and stellar population data suggest that the star 
formation history we've assumed is basically correct. 



While the stellar population models are largely insensi- 
tive to how stars formed early on, if there is a large burst 
at the end of star formation, it would lead to artificially 
young estimates if modeled with a quenched star forma- 
tion history. In fact, this can be seen for NGC 4522, 
where 2% of the stellar mass has to have formed in a 
burst at the end of star formation in order to match the 
stellar population data. The overall agreement of ages 
(even including NGC 4522) indicates that these galaxies 
only experience, at most, modest bursts at the end of 
their star formation epoch. 

Finally, the agreement of our quenching time and the 
timescales predicted by gas simulations imply that star 
formation is cut off when the HI is stripped. The 
value of tq tells us how long it's been since star formation 
has ended at the location where we measure the stellar 
populations. In the case of our sample of Ha truncated 
galaxies, we have measured tq just beyond the spatial 
truncation radius. Because of its relative proximity to 
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Fig. 12.— Same as Figure 14, but for NGC 4569. 
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the galaxy center, the region just beyond the truncation 
radius was the most difhcult to strip of any region already 
stripped by the ICM interaction. This implies that this 
region was stripped by the strongest pressure the galaxy 
had yet experienced. In the case of galaxies that are past 
peak pressure (as we believe to be the case for all of our 
sample), the region just beyond the truncation radius 
was stripped when the galaxy was experiencing peak or 
near-peak ICM pressure. The agreement between the 
time elapsed since peak pressure (from the ram pressure 
simulations) and the quenching time of star formation 
(iq from our work) means that star formation stops near 
the time of peak pressure. This agreement implies that 
gas stripping happens quickly and the stripping of the 
neutral gas (HI) is roughly coincident with the end of 
star formation. 

6.3. Locations of Stripping in Virgo 



In a study of nearby clusters, Solanes et al. (2001) find 
that, while HI deficiency extends out to 2Ra, the pro- 
portion of gas-poor spirals increases toward the center of 
the clusters. Additionally, Giovanelli & Haynes (1985) 
find that the HI deficient galactic fraction correlates with 
cluster X-ray luminosity in the sense that those clusters 
with higher X-ray luminosities (and presumably higher 
mass) have more Hl-deficient galaxies. 

It appears, therefore, that cluster processes are driv- 
ing HI deficiency and pushing gas out of galaxies. If it is 
the case that galaxies are only (or primarily) stripped in 
the core, this implies that there should be a relationship 
between galaxies' location in the cluster (Figure 18) and 
time since star formation quenching. While projection 
effects limit our understanding of the true distance be- 
tween a galaxy and M87 (at the core of Virgo), we can 
set lower limits to the clustercentric distance by using 
the projected distance from M87 (dg?). These are only 
lower limits because the galaxies are also displaced from 
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Fig. 13. — Same as Figure 14, but for NGC 4402. Note that H(5 is not detected, so there is no index diagnostic diagram shown. 



M87 by an unknown amount along the line of sight. 

In Figure 19, we show the distance of the galaxies from 
M87 plotted against tq for each galaxy. The dashed line 
in Figure 19 shows the distance a galaxy could travel in 
a given time if its velocity was 1500 km/s in the plane of 
the sky. The average line-of-sight velocity of a galaxy in 
Virgo relative to the cluster is ~ 700km/s, so 1500 km/s 
is a generous estimate for most galaxies. Any galaxies 
to the right of the line in the figure are either too far 
out to have been stripped in the cluster core or have un- 
usually large velocities in the cluster. This is simply a 
matter of travel time. Assuming some reasonable veloc- 
ity, it is impossible for the galaxies to the right of the line 
to have reached their current location in the time since 
star formation ended. Therefore, they can not have been 
stripped in the core. While several galaxies are consistent 
with being stripped in the core, there are three galaxies 
clearly to the right of the line and one galaxy marginally 
to the right of the line; we will take them case by case. 



6.3.1. NGC 4522 

NGC 4522, despite being far from the cluster core, has 
an exceedingly young population that precludes it from 
being stripped in the core. In addition to its young popu- 
lation, the large amounts of extraplanar HI and Ha, and 
enhanced polarized radio continuum emission and flat ra- 
dio spectral index along the leading edge all point to on- 
going strong ram pressure. Simple estimates of ram pres- 
sure are too low by a factor of 10 to adequately explain 
the observed stripping (Kenney et al. 2004). However, 
the "simple" estimates of ram pressure assume a smooth, 
static ICM. The existence of galaxies like NGC 4522 sug- 
gest that the ICM is either dynamic or "lumpy." Ob- 
served ICM shocks (Shibata et al. 2001) suggest that the 
ICM may be moving with significant velocity in the re- 
gion of NGC 4522. In Growl & Kenney (2006), we include 
a more complete discussion of the stellar populations of 
this stripped spiral galaxy. 
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Fig. 14.— Same as Figure 6, but for NGC 4064. 
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6.3.2. NGC 4580 

NGC 4580 is quite far from the cluster core but has an 
intermediate-aged outer stehar disk (iq — 475 Myr). In 
many ways, NGC 4580 is a prototypical stripped spiral, 
with normal star formation within a well-defined radius 
and no star formation beyond, but with an undisturbed 
stellar disk. One slightly peculiar aspect of this galaxy is 
the presence of stellar spiral arms beyond the truncation 
radius. These stellar spiral arms in the outer disk are 
stronger than in any of the other sample galaxies with 
comparable quenching times. Without gas, the spiral 
arms should dissipate relatively quickly, perhaps within 
a few galaxy rotation periods, but it appears these arms 
have persisted. 

6.3.3. NGC 4064 

NGC 4064 is very far out in the cluster, the furthest 
from M87 of any galaxy in our sample. Despite this, 



its stellar population points to star formation quenching 
only 425 Myr ago. It is possible that this galaxy was 
stripped of its star forming gas by an interaction with 
the ICM, but if that is the case, it did not happen in the 
core of Virgo. Optical imaging, kinematics, and CO in- 
terferometry (Cortes et al. 2006) suggest that this galaxy 
has recently undergone a minor gravitational interaction. 
However, this by itself is likely to not be enough to re- 
move the large amount of gas that this galaxy has appar- 
ently lost. The HI is truncated well within the stellar disk 
of NGC 4064, demonstrating that the gas was removed 
without significantly disrupting the stars. 

6.3.4. NGC 4424 

NGC 4424 is modestly below the 1500 km/s line and 
has a relatively young stripping age. This, coupled with 
the large HI tail observed in this galaxy (Chung et al. 
2007), suggests that this galaxy recently underwent an 
interaction. Optical morphology and stellar kinematics 
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(Cortes et al. 2006) demonstrate that NGC 4424 has re- 
cently experienced a significant gravitational interaction. 
Such interactions are typically more common in the clus- 
ter outskirts, where the relative velocities of galaxies are 
lower. However, the large HI tail pointing toward M49 
(Chung et al. 2007) suggests that a collision with that 
galaxy may have played a role in the stripping of the 
HI from NGC 4424. While this is clearly not the only 
environmental effect this galaxy has experienced, it ap- 
pears that a ram pressure interaction did play a role in 
quenching the star formation in NGC 4424. 

6.3.5. Galaxies Stripped in the Core 

The remaining galaxies fall near the vt — 1500 km/s 
line. This implies that, within the regions probed by 
our spectroscopic observations, star formation was in- 
terrupted on a recent core passage. If the galaxies had 
been stripped on a previous passage through the Virgo 
core, they would fall to the left of the line shown in the 



figure; that is, they would be closer to the cluster core 
than their stripping age and a travel speed of 1500 km/s 
would imply. The observation that these galaxies are at 
the minimum age that their distances allow suggests that 
the current appearance of these galaxies, with truncated 
HI and Ha disks, may be a short-lived state. One pos- 
sibility is that these galaxies are stripped further on a 
subsequent core passage. Note that because tq has only 
been determined just beyond the truncation radius, we 
cannot comment on how much of the outer disk was re- 
cently stripped; it's possible that only a narrow annulus 
of gas was stripped tq ago. As a galaxy enters a cluster, 
much of its gas is likely stripped on its first core passage 
but models suggest that, in a moderate-sized cluster like 
Virgo, it's unlikely that all of the gas will be stripped 
from an 0.5L* or l.OL* galaxy on that first passage (e.g. 
VoUmer et al. 2001). However, on subsequent orbits, gas 
at smaller and smaller radii may be stripped. This is 
particularly likely if cluster conditions and galaxy orbits 
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change, which are both Hkely in an unrelaxed cluster. 
However, our resuhs appear to show that some fraction 
of star-forming gas was truncated on these galaxies' last 
passage through the cluster core. 

6.4. Are Stripped Spirals "Failed" Passive Spirals? 

With the exception of small star- forming disks in their 
centers, the observed stripped spirals in our sample are 
similar to the "passive spirals" observed at higher red- 




FlG. 18. — The positions of the sample galaxies compared to 
locations of hard X-Ray emission (0.4-2.4 keV), as traced by images 
from the ROSAT X-Ray map of Virgo Cluster (Bohringer et al. 
1994). Two sample galaxies are not shown here; NGC 4064 is off 
of the figure to the north-west and NGC 4580 is off of the figure 
to the south. 
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shift (Dressier et al. 1999). It may be that Virgo is not 
massive enough to entirely strip the ISM from the spi- 
ral galaxies we have observed. However, in the largest 
clusters, the ICM density is higher by a factor of ~ 10 
and the relative velocities are higher by a factor of '-^ 2; 
this means that ram pressure is higher by a factor of 
r-^ 40 in these clusters. Therefore, it may be that the 
true passive spirals are only found in higher mass clus- 
ters such as those of Dressier et al. (1999) and that only 
truncated spirals are found in more modest clusters like 
Virgo. While our sample is not a complete sample of 
truncated spirals in Virgo, it is still interesting that none 
of the stripped spirals in our sample were stripped more 
than 500 Myr ago. If this observed state (truncated 
star-forming disk) were long-lived, we might expect to 
see stellar populations with a range of ages; from very 
newly stripped disks (i.e. NGC 4522 or NGC 4388) to 
disks that had been stripped several gigayears ago. Those 
galaxies with the most regular, symmetric Ha disks tend 
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to have the oldest values of tq (i.e. NGC 4405, IC 3392, 
NGC 4580), but those values are all less than or equal to 
500 Myr, a relatively short time compared to the orbital 
time in the cluster. We do not know whether the star 
formation in these galaxies will fade or the star forming 
gas will be completely stripped on a subsequent core pas- 
sage, but from the stellar populations of the outer disk, 
it appears that their current state is not long-lived. 

6.5. Metallicities of Galaxy Populations 

It is generally understood that galaxies create metals 
through stellar evolution processes. This means that, 
broadly, there is an age-metallicity relationship: in a 
given galaxy, the stars that form later tend to be more 
metal rich than the stars that formed first. While this 
is broadly true, the properties of any one region of a 
galaxy can be largely determined by the local star for- 
mation history. In order to avoid adding free parameters 
to the SB99 models, we summed stellar models at a con- 
stant metallicity; all stars in a given model (regardless 
of the age) have a fixed metallicity. This is, of course, 
a gross over-simplification of the stellar population of 
the galaxy, but due to the relative insensitivity of the 
Balmer indices to metallicity (particularly at the young 
ages that are most relevant for our study), we have cho- 
sen not to model any sort of enrichment history. As a 
result, the metallicities implied by the [MgFe]' index are 
the luminosity-weighted average metallicity for the entire 
galaxy and not the metallicity of the current generation 
of stars. It should be noted that the luminosity-weighted 
average metallicity is slightly sub-solar for most of our 
sample of spiral galaxies. This is expected, as galaxies 
with solar- metallicity current generation stars likely have 
a luminosity-weighted mean luminosity that is slightly 
sub-solar; in galaxies with realistic enrichment histories, 
the older stellar generations have fewer metals than the 
current stellar generation. 

7. SUMMARY 

We have presented observations of several gas-stripped 
spiral galaxies in the Virgo cluster. These are galaxies 
that have apparently been partially stripped of their gas 
by interactions with the ICM, but all still have ongo- 
ing star formation in their centers. We present optical 
spectroscopy that demonstrates that the outer disks were 
stripped of their star-forming gas within the last 500 
Myr. We find that, while some galaxies are consistent 
with being stripped in the cluster core, others appear to 
have been stripped outside the core. There is evidence 



that some outer cluster galaxies have undergone gravi- 
tational interactions, but others appear consistent with 
ram pressure stripping. These galaxies provide evidence 
that the ICM is not static and smooth. In at least one 
case (NGC 4522), there is independent evidence that the 
ICM pressure is stronger than simple estimates would 
predict. If the effects of ram pressure are important well 
outside the cluster core, the HI deficiencies observed at 
large radii in clusters (Solanes et al. 2001) can be ex- 
plained by ISM-ICM interactions. These interpretations 
imply that the "reach" of the ICM extends well beyond 
the cluster core and that ISM-ICM interactions may have 
a greater impact on morphologically-driven galaxy evo- 
lution than simple ICM models would suggest. 

Furthermore, for three galaxies that have been mod- 
eled by gas dynamics simulations, we find agreement be- 
tween our determined quenching times and those sim- 
ulations. These results imply that our galaxies experi- 
enced essentially constant star formation prior to being 
quenched and that there was, at most, a modest star- 
burst at the time of quenching. It also suggests that star 
formation is effectively halted by these strong ICM inter- 
actions and molecular gas is effectively stripped quickly 
by ICM pressure. 

Taken together, our observations show that gas strip- 
ping is an important process in environmentally-driven 
galaxy evolution. As ISM-ICM interactions are capa- 
ble of rapidly halting star formation, it is reasonable 
to suppose that a significant passive galaxy population 
may evolve from spiral galaxies in clusters. While there 
may be multiple paths to SO galaxies, it appears that, 
in Virgo, we are observing a process that can efficiently 
halt star formation in cluster spiral galaxies. 
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